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Non-minimally coupled theories of special potentials are analyzed numerically. Such theories yield 
equations of state uo < — 1 and oscillations of the cosmological expansion, which are favored by the 
recent analysis of observations. Fitting these theories to the Gold Snla dataset, we obtain results 
comparable with other models. A potential of the form V((j>) = Voe ai ^ yields Xmin = 170.127. 
Similar results are obtained for potentials of the form V(<f>) — Vo + ai<j> n . 



I. INTRODUCTION 



It has been suggested strongly by detailed observations of distant Type la Supernovae (Snla) that our universe 
may be in a phase of accelerating expansion. On the other hand, measurements of cosmic microwave background 
(CMB) 0] and other surveys Q indicate the spatial flatness of the universe. The simplest explanation for these obser- 
vations is a dominating cosmological constant, whose equation of state is lo\ — p\/ p\ = —1, with a complementary 
Cold Dark Matter (LCDM): 

H 2 (z) = H$[flA + n m0 (l + z) 3 } (1.1) 

where z is the redshift, Q m o is the present value of matters (including the dark matter, whose existence has well been 
confirmed) and Q\ that of the cosmological constant. Spatial flatness demands that J7a + f^mo = 1. 

LCDM is a simple proposition but requires extreme fine tuning of the cosmological constant p\. More so, its fit with 
the Gold dataset (157 most updated and reliable set of Snla) is not very good. Alternative models with evolving dark 
energy, such as quintessence, phantom, extended gravity, were then put forward Q, most of which contain a rolling 
scalar field. On the other hand, recent analysis of Gold dataset and other observations seem to favor an equation of 
state u < — 1 and oscillations of cosmological expansion 0, |(| . 

In normal theories of scalar fields, the equation of state is usually 

_ ^ - m 

u » - wrwy (L2) 

Obviously, it is not possible to have ui < — 1. One attractive way out of this situation is to consider extended theories 
in which scalars are non-minimally coupled. These extended theories appear naturally and widely in string theories 
and theories of extra dimensions in general. 

In this paper, we consider simple theories of a non-minimally coupled scalar field <p which is oscillating around the 
minimal of the potential V(^). In this type of theories, u) < — 1 and oscillations of the cosmological expansion can 
be obtained readily, as the oscillating scalar field is now strongly coupled to the gravity, and it is well known that 
non-minimally coupled theories can violate energy conditions. 

Fitting to the Gold dataset, these models yield sensible results comparable with other models. The best fit param- 
eters of these models show that the oscillating field <f> may account for both the dark matter and the dark energy in 
the cosmological expansion. The behavior of <f> depends on the energy density, so </> could differ significantly within 
or without galaxies. Properly averaged over space, <fr is scale-dependent. We have concentrated our analysis mainly 
on large scales, though it is possible that the dark matter appearing on other observations is also due to motion of 
the same scalar field. Interestingly, our models which have properties of lu < — 1 and oscillating cosmological expan- 
sion seem to have better fit, which can be interpreted as circumstantial evidence for oscillations of the cosmological 
expansion 0. For comparison, ll H 2 (z) — z" and i( u>(z) — z" relations are plotted for our models, LCDM model and 
for the best fit model OA Var.(l) in Q. 

The paper is organized as follows. In section II we present fundamentals of non-minimally coupled theories and 
outline the procedure of fitting models to the Gold dataset. We then analyze two specific models for illustrations, 
with results presented in detail in section III. In section IV, we conclude with some discussions. 
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II. FUNDAMENTALS OF NON-MINIMALLY COUPLED THEORY AND THE FITTING PROCEDURE 



Generically, the lagrangian density for a non-minimally coupled scalar field theory assumes the following general 
form 



S = d^x^J—g 



\f(4>)R - ^W(c/>)d^d^ - U(<l>) + L N 



(2.1) 



where Ln is the total lagrangian density of fields and matters in the universe other than the field (/). In this paper, 
the theory is parameterized such that W(4>) = 1 and F((f) of the simple form 



F(4>) = 1 - £0 2 



(2.2) 



here we have set 8irG = 1. Assuming a flat Friedmann-Robertson- Walker metric, i.e., we assume a flat prior in all 
calculations of this paper, the Einstein equations are |7j,|8|]: 
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{ Pn +pn) + 4> 2 + f-h-f 



(2.3) 
(2.4) 



where /?jv and come from the Ljv, representing, respectively, the energy density and the pressure of matters and 
fields in the universe other than <j). In what follows, we will approximate these matters and fields as a perfect fluid 
with = 0. 

For simplification, one makes the following redefinitions 

Q = H 2 /Hl V{4>) = U(<j>)/H 2 (2.5) 
where H is the present value of expansion rate H. From Eqs (|2.3|> and l|2.4|) . one obtains p^Hof 
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F-F'{l+z)-^-{l + zf 



(2.6) 
(2.7) 

by those with 
(|23|l yields a 



where ft no — Pno/3Hq. In these two equations, we have replaced derivatives with respect to time 
respect to the redshift z, which are in turn denoted by primes. Substitution of Eq 12.7fl into Eq 
second-order differential equation which includes only field <j>. 

To be consistent with measurements in the solar system |ll| , stringent constraints have to be put on these theories 
Defining a post-Newtonian parameter 



7- 



(dF/d(/)) 2 



F + (dF/d(/)) 2 

solar system tests give the following (present) upper limit [Tll | 

| 7 - 1| < 4 x 1CT 4 , 

which is equivalent to 



1 

F 



dF 



< 4 x 10" 



|^(0)| < io- 



To generate suitable oscillations in the cosmological expansion, f should be of 0(1) (see the next 
constraint is now converted into a constraint on <fi(0). Numerically, (j>(0) < 10~ 3 is strict enough. 
Tests in the solar system constrain only the behaviors of <f> in our galaxy. According to its equation 



(2.8) 
(2.9) 

(2.10) 

section) . The 
of motion 







(2.11) 
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(where R is the Ricci scalar) the evolution of <j> depends on the energy density, for R is proportional to the energy 
density. Usually the potential V(4>) is tiny, such as the ones appearing in the next section. Within galaxies, where 
the density of matters is much larger compared with the average density of the universe, the potential term in the 
above equation can be neglected. Outside of galaxies or averaged over large scales, V(<p) plays the main role and the 
frequency of oscillations of <fi would be much smaller. So the magnitude and the evolution of <f> are scale-dependent 
and sensitive to positions, which produce varied effects at small and large scales. As we are interested in cosmological 
expansion at large scales, </>(0) is much less constrained. 

For simplicity, we will take </>(0) = as one initial condition. Other proper choices of 4>(0) will yield qualitatively 
similar results. On the other hand, the effective Newton constant 



G„fF — 



1 2F + 4(dF /d<j>) 2 



F2F + 2,{dF/d(j)) 2 
changes in time, which is also constrained by present observations 



Gcff 



< 6 x 10~ 12 yr -1 



(2.12) 



(2.13) 



2 = 



By choosing the particular initial condition 0(0) = 0, one automatically has 



Gcff 



G ff 



= 



2 = 



The second initial condition is obtained by evaluating Eq (|2.7I) at z = 0: 

0' 2 (O) =6 - 2V(0) - 6Q m 
From Eqs l|2.6|l . (|2.7|l and these initial conditions, one obtains the current equation of state 

.(o) = i^Ao)-M 



(2.14) 



(2.15) 



(2.16) 



If £ is large enough, oj(0) can be less than —1 easily. 

Given these two initial conditions, the second-order differential equation of 4> can be solved, at least numerically, 
once the exact form of V(<f) and the value of f2jvo are also known. Substituting the solution <p(z) back into Eq (|2.7|l . 
we get H(z) for arbitrary redshift z. Thus, we solve the problem completely. It is then straightforward to compare 
the theory with observations. 

To be definitive, one calculates the goodness of fit of the theory with the observed Gold dataset: 
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m° b (z k )~m th (z k ,H,M)y 



fe=i 



m° b (z k ) 



where 



is the apparent magnitude and 



m m (z, H,M)=M + 5 log 10 (1 + z) 



M = M + 51og 10 ( 



dz 



H(z') 



Mpc ' 



25 



(2.17) 



(2.18) 



(2.19) 



is the magnitude zero point offset. 

To find the best fit parameters of the model, one minimizes x 2 . Since \ 2 ls a second-order polynomial of M, the 
minimization procedure with respect to this parameter is straightforward and simple. The final result of Xmin can De 
expressed in terms of quantities independent of M : 



X 2 {H) = A-B 2 /C 



(2.20) 
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TABLE I: Best fit parameters and error bars for two types of potentials in non-minimally coupled (NMC) 
theories, compared with the LCDM model and the OA Var. model in [6]. The first two rows are for the first 
form of potential. Shown in the first row are the results with Qno as a freely fitting parameter, while shown 
in the second row are those by setting Qno = 0.045. The next three rows are for the second form of potential, 
with power index n — 2, 4, 6, respectively. When estimating the error for each parameter, all other parameters 
are fixed at their best fit values, except that f2jvo is fixed at 0.045 in the second and the fifth rows. 



Model 




V(cj>) or H 2 (z) 


V 


a i 






Xmin 


NMC E(0) 




V{4>) = V Q e a ^ 2 


2.18 ±0.08 


25.91^ 


1.94 ±0.67 


U.UOQ_q Q21 


170.127 


NMC E(l) 




V{4>) = V e a ^ 


2.19 ±0.08 


25.0±*;S 


i-79l° ;S 


0.045 


170.135 


NMC M2(2) 




V{<t>) = Vb + a 1( t> 2 


2.05 ±0.015 


-34.71*;? 


16.01?;° 


0.300ir o °ol 


171.870 


NMC M4(3) 




V(4>) = V + ai^ 


2.30 ± 0.08 


i446t2:S 


0.50 ±0.18 


0.078lg : °» 


172.494 


NMC M6(4) 




V{<j>) =V + aic/) 6 


n R7 +0.03 
z - D '-0.04 


i8.6±g:3 


0±0.04 


0.045 


173.963 


OA Var. (5) 


H 2 (z 
(1 + *) 3 | 


) = Hi [0.3(1 + zf + 0.7 + 0.13 
cos (6.83z + 4.57tt) - cos (4.577r)]] 










171.733 


LCDM(6) 


H 2 . 


[z) = [0.31(1 + zf +0.69] 










177.072 



where 



_ ™ (m ob (z k )-m th (z k ,H,M = 0)) 2 

k=l rn° b {z k ) 

¥1 (m ob (z k )-m th (z k ,H,M = 0)) 
B = 2 ( 2 - 21 ) 



C = 



k=l a m° b (z k ) 
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/ . 2 

k=l m° b (z k ) 



Eq l|2.20|l is the starting point of our numerical analysis. By minimizing this quantity, one obtains the best fit values 
of relevant parameters. As we shall see in the next section, both lo(0) < — 1 and oscillations of cosmological expansion 
are favored by experiments. 



III. NUMERICAL RESULTS 



Taking two special potential forms for illustration, we now calculate the fit of non-minimally coupled theory with 
the Gold dataset and find the best fit parameters in the process. In order to obtain acceleration and oscillations of 
the cosmological expansion simultaneously, we consider potentials that have non-zero minima Vq at the point = 0. 
Though we have only analyzed two forms of potentials in this paper, our analysis can easily be extended to other 
potentials. Other potentials of this type could yield comparable fitting results, as in general one can easily obtain 
appropriate oscillations and acceleration of the cosmological expansion from this type of models. All results are shown 
in Table I. The errors are obtained with an increase of Xmm by 1, while we have assumed that the value of f2jvo is 
at lest 0.045, which is the density of all know matters including baryons and CMBR. When estimating the error for 
each parameter, all other parameters are fixed at their best fit values. For comparison and concreteness, we will plot 
the evolving of <f>(z), H 2 (z) and u>(z) for all models in FIG. [SI FIG. El and FIG.0 respectively. 
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A. V{<t>) = Voe^ 2 

This potential appears naturally in supergravity theories (SUGRA) |l3|. This form was also reconstructed from 
Snla Gold dataset in , where it is only a good fit of the polynomial parametrization model in reference [(| and <f> is 
not oscillating. In total, we have four free parameters: £, &no, Vq & n d (X\. Their best fit values are shown in first line 
of Table I (NMC E(0)). 

FIG. n shows the dependence of Xmin on ^-nq, when other parameters are allowed to vary freely. The best fit of 
£Ino is about 0.066. However, Fig. ^ indicates that one gets almost no change in % 2 if Qno is set to be 0.045. The 
latter is the value of all known matters including baryons and CMBR [l4|. That is to say, the existence of a single 
non-minimally coupled scalar field <f> may play both the roles which dark matter and dark energy would play in LCDM 
on the cosmological expansion. Of course we can split this effect into two parts: dark matter (which also generates 
oscillations of the cosmological expansion) and dark energy. Clearly, this splitting is rather artificial, as there are no 
obvious distinctions between these two parts in this model. They can just be regarded as a whole as the expansion 
of the universe is concerned. Now that there is no isolated dark energy we plot in FIG. the total equation of state 
of the universe. Dark matters from other observations can also come from this field [l^, though the behavior of the 
field should be scale-dependent ^(| . In this paper we only consider the effect of the field 4> on cosmological expansion 
at large scales. 

Clearly the bigger Q,no is, the worse is the fit to Gold dataset. This may be seen as follows. To obtain enough 
amplitude of oscillations of field <f>, i.e., to obtain enough oscillations of cosmological expansion to fit observations, 
</>' 2 (0) should be large enough. On the other hand, to have a dominating dark energy, Vq should be positive and large 
enough. According to Eq i|2.15[l . fijvo cannot be too larger. In the rest of this subsection, we will fix ft no = 0.045 to 
analyze other parameters. A re-analysis of this model with this value of fijvo is given in the second row of Table I. 
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FIG. 1: The dependence of x 2 on Qno with other parameters freely, and the dash line is for Xmin + 1- 

The parameter Vq here plays somewhat the role of cosmological constant, with the effective 51a = Vq/3. Our best 
value Vq ~ 2.19 is consistent with the best fit value JIa = 0.69 in LCDM model. If Vq departs too much from the 
best fit value, we will have a comparatively large x 2 ■> which is confirmed from FIG.|3 In FIG.|2]we have plotted the 
X 2 — Vq relation by fixing the other parameters as: SIno — 0.045, £ = 1.79 and a\ — 25.0. 




v 

FIG. 2: x 2 v - s - Vo with other parameters fixed as: Qno = 0.045, a\ — 25.0, £ = 1.79, and the dash line is for Xmin + !• 
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Now we turn to the parameter £, which represents the coupling strength between gravity and <j>. As </> oscillates 
around the minimum of potential V(cj>), £ affects the amplitude of oscillations of cosmological expansion. If £ is very 
small, there will be almost no oscillation appearing in the expansion rate H 2 and it is difficult to cross the u> = — 1 
line. Its behaves nearly the same as that of the minimally coupled theory @|. If £ is very larger, the motion of the 
universe will be drastically changed. So the value of £ should be moderate. In fact, the best fit value of £ is of the 
order 0(1) (sec Fig. ©). 
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FIG. 3: x 2 v.s. £ with D. N o = 0.045, Vo = 2.19, di = 25.0. The dash line is for xLi„ + 1- 

Last but not the least, the parameter a\ can affect frequency of oscillations, its effects on \ 2 are shown in FIG.0] 
As mentioned above, we have plotted 4>{z), H 2 (z), and u>(z) in terms of z for all models, in Figs. OIBJandEI One 
clearly sees oscillations of <f> and of the cosmological expansion over time in this type of model. The universe may well 
be just in a period when the lu is in the valley and later on will transit to a period of deceleration slowly. 
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FIG. 4: x v.s. oi with Cl N0 = 0.045, Vo = 2.19, £ = 1.79. The dash line is for xLin + 1- 



B. V{4>) = V + ai<j> n 

This form of potentials can be regarded as an extension of a constant dark energy, i.e., a cosmological constant Vo 
is added with a field <f> of power potential. The non-minimally coupled scalar field mimics the effects of dark matter 
and produces oscillations in the expansion rate H 2 . We have considered cases of n = 2, 4, 6, with their best fit values 
of parameters shown in Table I. As shown in FIG. ©and FIG. the expansion patterns of these models are similar 
to those of LCDM and oscillate around the line of LCDM. 

For n = 2, the best fit value a\ < and the potential a\4> 2 is negative, but coupling of the scalar field with gravity 
can drive up the potential. The field (f> is still oscillating around <f> — if all parameters arc within certain range 
around the best fit values. In this model the constraints on parameters arc quite strict, as shown in Table I. In Figs. 
03 E]and0 we see that oscillations of the cosmological expansion only appear near the cf> = points, while in other 
places there is essentially no oscillation. 

For n = 4 and more so for n = 6, the cosmological constant Vo is quite large. Then there is not much dark matter, 
which seems to be in conflict with other observations of dark matter. It is possible that we can have enough dark 
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matters at relatively small scales, as the effect of <j) is scale-dependent. For n = 6, the best fit value of £ is 0, but there 
are still oscillations in the expansion rate H 2 , with a much smaller amplitude. Its best fit value of flj^o is negative, 
much smaller than 0.045, the value of known matters. We set the value Qno = 0.045 in Table I by hand. 

IV. CONCLUSION 

We have numerically analyzed the expansion of the universe in non-minimally coupled scalar field theories with two 
special forms of potentials: Voe ai ^ and Vb + ai<p n . These types of potentials have a no-zero minimum at <j> — 0, which 
generates the necessary acceleration of the cosmological expansion. On the other hand, the oscillating <f) produces 
oscillations of the expansion and yield an equation of state w < —1, as favored by the fit with 157 observed Snla 
Gold dataset. We have further compared our models with the Gold dataset to find the best fit parameters. The best 
fit potential is of the form Voe ai ^ with Xmin = 170.127. The best value of all matter density in universe except (/) 
is found to be fi^o — 0.066, which is very close to all energy density of the known matters including baryons and 
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FIG. 5: Field <j>(z) v.s z in NMC models. 




NMC M2 ( 2 ) 



FIG. 6: Expansion rate H 2 {z) v.s. z in six representative models. Here H(0) has been normalized to be unity. 
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CMBR. So the non-minimally coupled scalar field <f> may provide the main source of dark matter as well as dark 
energy. Of course further investigation will be needed to make these assertions more definitive. 
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